A POSSIBLE EXPLANATION WHY tb± ~ r^o BUT r^± ~ 2t^o 



Hui-Shi Dong^, Xin-Hcng Guo"^'^, Xue-Qian Li^'^, and Rui Zhang^ 



1. CCAST (World Laboratory) P.O.Box 8730 Beijing 100080, China. 

2. Department of Physics, Nankai University, Tianjin, 300071, China. 

3. Department of Physics and Mathematical Physics and 
Special Research Center for Subatomic Structure of Matter, University of Adelaide, SA 5005, Australia. 

4. Institute of High Energy Physics, Academia Sinica, Beijing 100039, China. 



Abstract 



Data show that Tg± ^ r^o, but r^ji ^ 2r£)0. The naive interpretation which 
attributes td± ~ Sr^io to a destructive interference between two quark dia- 
grams for decays, definitely fails in the B-case. We investigate Close and 
Lipkin's suggestion that the phases for producing radially excited states ^p2s in 
the decay products of B-mesons can possess an opposite sign to the integrals for 
ipis decay products. Their contributions can partially compensate each other 
to result in Tg± ^ Tga . Since D-mesons are much lighter than B-mesons, such 
possibilities do not exist in D-decays. 



PACS numbers: 13.20.Fc, 13.20.He, 12.39.-x, 13.25.-k 



A Possible Explanation Why r^i ~ r^o But r£)± ~ 2x^)0 



I. Introduction. 

The naive explanation for r^ji ~ 2t£)0 is that a destructive interference between two quark 
diagrams for reduces the strength of decay amphtudes and thereby elongates the life of 

. More explicitly, if the lifetime of a meson is mainly determined by the Cabibbo favored decay 
modes, for Z?+ there is only one topology £)+ K'^M^, (where M generically refers to 7r,p etc 
and K to strange mesons) whereas for there are two channels 1?° K^M^ and M^. 
For the Z?+ decays, the two quark diagrams shown in Fig.l (a) and (b) interfere, while for the 
two diagrams (c) and (d) correspond to two different modes, so do not interfere. For the B decays, 
similar diagrams exist and there could be also destructive interference in B~ decays. However, the 
experimental data show that tb± ^ Tgo 

The explanation for the lifetime differences in D and B cases invloves nonperturbative QCD 
phenomena. Actually some authors Q| proposed the so-called Pauli Interfernce (PI) mechanism 
as a correction to "pure" spectator mechanism for taking into account of the light degrees of 
freedom. The PI effects only exist in and decays but not in and decays. Based 
on QCD, Bigi et al. introduced a virtual gluon so that one of the quark produced by the 
weak decay of the heavy quark interferes with the spectator quark. In this mechanism, the PI 
term modifies the "pure" spectator diagram and it is found that such interference is destructive 
and is proportional to Fo/mg (Q=b or c). This mechanism partly explains why r^ji ^ 2rpo and 

In the present work we try to investigate the lifetime differences in another way which is based 
on the idea of Close and Lipkin. Recently Close and Lipkin[|j have analysed the data on low 
lying exclusive quasi-two body final states in both D and B decays. They noted that in D decays 
the sign of interference in exclusive channels is still amiguous while in B decays there is a clear 
and uniform tendency towards constructive interference between the color favoured and colour 
suppressed exclusive channels where all final state mesons have nodeless wave functions. They 
noted that in B decays, in order that r^i ~ r^o H, this interference must be compensated in as 
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yet unmeasured channels. They suggested that the sign of interference may be changed in channels 

where excited states of the decay products, whose wavefucntions contain nodes, are involved. It 

is the motivation of the present paper to include the contributions from the excited states of the 

B decay products so that constructive interference is obtained in decays. Such excited states 

only excist in B decays but not in D decays because of phase space requirement. It will be shown 

that in our model the lifetime differences in B and D mesons can also be explained. 

The effective Hamiltonian of non-leptonic decays in the D-case is 

Gf - - 

Heff = -j=V*^Vud[cis-i^,{\ - 75)cM7''(l - 75)d + C2S7^(1 - 75)dw7^(l - 75)0], (1) 

where ci — and C2 = '^^ 2^ ' ^^"^ renormalization group equation (RGE) we have 

_ ^ as(TOg) ^12/25^ as(mg) ^^2/23. ^ _ 1 



With the Fiertz transformation, the coefficients ci and C2 in eq.(jl|) should be replaced by oi and 
02 with 

fli = ci + ^C2, and 02 = C2 + ^ci, (3) 

where ^ is 1/Nc if the factorization assumption holds perfectly, otherwise ^ = (1 + 5)/Nc where S 
denotes a color-octet contribution proportional to < A^'A" > [^[^|^. Recently, Blok and Shifman 
gave a more theoretical estimation ||l0| , but they also pointed out that the obtained value is not 
accurate for practical calculations. Generally, ^ is a negative number ranged between to —1, so 
that ^ takes values between to l/N^- Later, in our numerical calculations we will take 6 as 0, 
—0.5 and —1 respectively. 

For the B-case, we have a similar Hamiltonian as eq.(^ 

H,ff = ^v,bV:,[cfhj^{l - 75)M7^(1 - 75)" + ci^'hj^il - 75)^7^(1 - I5)b] (4) 



and coefficients €{ and a 
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^(.) ^ (i^),2/23. ef = ^, (5) 

as[M^) ^ Lb 

whereas al ,03 have similar forms in analog to that for the charm case. 

It is noted that in the case of D meson decays ai is positive and and 02 is negative. From 
the data of D-physics the value of 02 is about —0.5 [O . In 13+ decays, ai term corresponds to 



the external W-emission while a2 to the internal W-emission, naturally a destructive interference 
would occur between the two quark diagrams. 

In the following wc will express the corresponding transition amplitudes as Ai and A2 which 
arc proportional to ai and a2 respectively, thus Ai = kiUi and A2 = K,2a2 where ki and K2 are 
the hadronic transition matrix elements. 

Then we have the amplitude square as, 

I < K\+\H,ff\D+ > |2 = + \A2\^ + 2i?e(AiA*). (6) 

Beside a common phase factor such as the Cabibbo-Kabayashi-Maskawa phase, both Ai and A2 are 
real. Thus if Ai ■ A2 is negative, this is a destructive interference. Otherwise we have constructive 
interference. 

In contrast, for D'^ decays, 

< K-'!i+\Heff\D^ >(xai and < K^T:°\Hef f\D^ >(x a2. 

We can roughly assume 

<K~Ti+\H,,ff\D° Ai and < i?"7r°|i7e//|£'" >~ ^2- 

Thus if we only consider the C-K-M favored channels which dominate the lifetime of D-mesons, 
we have 

T[D+) = + + 2i?e(^a*)) X (7) 

r(i?°) - {\A^\^ + \A2\^)xLIPS, (8) 

where LIPS is the Lorentz-invariant-phase-space of the final products. If A2 ^ — 0.26Ai, one can 
numerically obtain T{D'^) ^ 2T{D^) (or T(£)±) ^ ^tjjo). Of course, the other channels (Cabibbo- 
suppressed) and semi-leptonic decays all contribute to the lifetime, so this obtained number is not 
rigorous. However, since the Cabibbo favored channels dominate, one can expect that a solution 
for Ai and A2 does not deviate much from the aforementioned value. 

Taking as(M'^) = 0.118 [|l|, one can obtain a ratio of A1/A2 for D-decays to be roughly 
consistent with the required value. By our recent knowledge, the hadronic matrix elements can be 
evaluated more easily in terms of the heavy quark effective theory (HQET) fl^ j. 
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In the same scenario and by eqs.(^, ai and 02 are still of opposite signs in B-decays. It 
is a consequence of renormalization group equation which is proved to be valid for perturbative 
QCD. If so, one could expect a result similar to the D-casc that Tg- ^ 2tbo. However, this does 
not coincide with the data for B-decays. 

The lifetime is very close to that of B° as r^ci) (1.62 ± 0.06) x lO'^^ g ^^^^ ^ 

(1.56 ±0.06) X 10^^^ s [|l[. There could be small measurement uncertainty as Tj^i) ~ 1.47 x 10~^^ 



s, r(BO) - 1.25 X 10^12 by the ALEPH collaboration ||13|ll5| and T(s±) ~ 1.72 x lO^^^ 
T(so) - 1.63 X 10^12 by the DELPHI collaboration Hi. 

Similar quark diagrams exist in B-dccays, namely there are both external and internal W- 
emissions for D^tt^ which destructively interfere, but for D^tt^ and B'^ D^tt^ 

corresponding to external and internal W-emissions respectively do not interfere. Thus if that is 
the case, one would wonder why r^i is so close to r^o. 

To fit the data of B decays, one needs to take a positive value for 02 [0. This contradicts 
to the result of RGE which is obviously correct by the perturbative QCD theory and there is no 
doubt of application of perturbative QCD at the nib energy region. 

However, one can notice that even though Ai, A2 are proportional to ai, 02 respectively, they 
also possess certain factors corresponding to the hadronic matrix elements. These hadronic matrix 
elements involve some overlapping integrations of the decay parent and daughter wavefunctions. If 
the integrations can contribute a negative sign, the interference between two diagrams would turn 
over to be constructive and it may be equivalent to an "effective" positive 02 value. 

The hadronization process is very non-perturbative and we cannot evaluate it accurately, so 
that we attribute the non-perturbative effects into the parameters of meson wavefunctions which 
exist in the overlapping integration. To evaluate such overlapping integrations, one needs to invoke 
some concrete models and later we employ the non-relativistic quark model. Since the decaying 
B-meson is a pseudoscalar at -01^ radial ground state, if the decay product is at "tpis state, the 
overlapping integration would certainly be positive, however, if the decay products can be radially 
excited states V'2s, the integration can turn sign (see next section for details). Because D-meson 
is much lighter than B-mcson, it does not seem to exist ip2s states as decay products of D, but 
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definitely there should be 'ip2s excited states showing up as decay products of B-meson. This 
change may modify the whole picture and finally leads to a consequence that Tq{±) ^ Tgo. Later 
our numerical results will show that the involvement of the "023 decay products can indeed do the 
job. 

In next section, we give the formulation in every detail and in Sec. Ill, we present our numerical 
results while the last section is devoted to conclusion and discussion. 

II. Formulation 

(i) The transition amplitudes. 

As usual, we ignore the W-exchange and annihilation diagrams because the two fast quarks 
would pick up a quark-pair from vacuum and speed them up Q. Even though the factoriza- 
tion approach is not very reliable to evaluate the internal W-emission diagrams, we may use a 
phenomenological parameter 6 to compensate it. Therefore by the vacuum saturation 

< K-Tr+\ai{sc){ud) + a2{sd){uc)\D° >= 

ai < TT+\{ud)\0 >< K-\{sc)\D° > +02 < K-TT+\{sd)\0 >< 0\{uc)\D° >= 

aiUp^ < K-\{sc)\D'' > +a2.fDp'i < K-Tr+\{sd)\0 > (9) 

where {qq') = q^fii^ — 75)9'- The second term corresponds to an W-annihilation diagram and 
obviously is much smaller than the first one as it is proportional to foin^^ ~ '^^)- argued in 
literatures this term is negligible and we will omit such contributions in later calculations. Then 
we also have 

< X°7r°|ai(sc)(ud) + a2{sd){uc)\D° a2/KP^^ < tt°\{uc)\D° >, (10) 

and 

< K\+\ai{sc){ud) + a2{sd){uc)\D+ >^ a^Upi^ < K°\{sc)\D+ > +a2fKp'ii < Tr+\{uc)\D+ > . 

(11) 

Instead, for P ^ PV 

< K-p+\ai{sc){ud) + a2isd){uc)\D^ > = aifpUipe*'' < K-\{sc)\D° >, (12) 
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and 



< K-*T:+\ai{sc){ud) + a2isd){uc)\D" > = ai/^p^ < K-*\{sc)\D° >, (13) 

< K°*TT°\ai{sc){ud) + a2{sd){uc)\D" > = a2fK'mK*e*'' < 7r"|(uc)|L»" >, (14) 
< ifV|ai(sc)(wd) + a2{sd){uc)\D" > = a2fKp'ii < P°l(Sc)|i?° >, (15) 

< K°*TT+\ai{sc){ud) + a2isd){uc)\D+ > = aif^p^^ < K°*\{sc)\D+ > + 

a2fK'e*^,mK' < Tr+\iuc)\D+ >, (16) 

< K°p+\ai{sc){ud) + a2{sd)iuc)\D+ > = ajpe^'^rup < K°\{sc)\D+ > + 

a2fKp';i < P+\{uc)\D+ > . (17) 

For P^VV 

< K°*p°\ai{sc){ud) + a2{sd){uc)\D^ > = a2/x*mK.e*'' < p°\{uc)\D^ >, (18) 

< K-*p+\ai{sc){ud) + a2isd){uc)\D° > = aifpTJipe*" < K-*\{sc)\D° >, (19) 

< K^*p+\aiisc){ud) +a2{sd){uc)\D+ > = aifpei^nip < K°*\{sc)\D+ > + 

a2fK'e*j^,mK' < p+\{uc)\D+ > . (20) 

The above formulae indicate that the external and internal W-emissions in D'^ decays interfere. 
For the B-case, we can have similar expressions with an effective hamiltonian eq.(^) and corre- 
sponding coefficients af , af in eq.(^. 

(ii) The matrix elements. 

It is noted that in the scenario of factorization, the hadronic matrix elements are related to a 
weak transition ||l^ , for P ^ P 

f2 _ /i,,r2 



and 



< X\j,\I [Pj + Px- ^^l-^q),F,{q^) + ^^LJ^q^F,{q% (21) 
with q = Pi-Px and Pi(0) = Po(0). For P -> T/, we have 

< > = MiImx^ ""''*"^'^'''^^'^^^ + z^2Afx.9pAo(92) + 



e* • q 



g2 



2Mx.g^A3 (22) 



with ^3(0) = Ao(0) and here 



So our task is to calculate the form factors. Taking the nearest pole approximation 

^1(9') « Y^^T]^ for - Fx, (24) 
ng')^ , ^^..2 ^ ^o(g') = ^-%T72 for P/^Px., (25) 



where Afi, M2, M3 are masses of mesons corresponding to the nearest poles which can be found in 
the data book. With this approximation, to evaluate the form factors, one only needs to calculate 
the constant parameters Hq = hi, hy, hAn B^nd hA^ — hAo, which turn out to be the values of 
the form factors at the unphysical kinematic region = and we will use the non-relativistic quark 
model to calculate them. Moreover, for the case of a pseudoscalar B or D meson transiting to a 



vector meson, we use the helicity amplitude method |17 which can much simplify our calculations. 

The parameters are related to an overlapping integral over the wavefunctions of initial pseu- 
doscalar and final pseudoscalar or vector mesons. To carry out the integration, one needs to invoke 
concrete models and the most popular one is to take the wavefunction of harmonic oscillation 
potential as the orbital part. In the Bauer-Stech-Wirbel approach the following wavefunction 
model is employed 

-2 1 ml ~ ml 



R^{pT,x) = iV„y^(T^exp(~p2,/2^2) .g^p(__(^_ _^^__^)2) (26) 

where Nm is the normalization factor while Guo and Huang [ll6| used the following wavefunction 
form in the light-cone formalism 

Rra{x,k,,rt) = Ae^^bH^^^^ + ^ll^)) . (27) 

These wavefunctions apply in the infinite-momentum frame. Here instead, we choose the wave- 
function at the rest frame of the decaying meson |T^. Everything in the picture is non-relativistic, 
but it is accurate enough for the qualitative conclusion and we will discuss it in the final section. 



Here we only list the radial wavefunctions of tpis and V'2s and the others can be found in ref. p9| 
i;u - (^)i/2exp(-i/32r2)V^roo(0,0), (28) 

y/TT 2 

where adding a factor ^/m is for proper normalization, and 

i,,, = {-^Y'^ii - 2/3^2) exp(--/?2r2)roo(0, <^)V^, (29) 

where [3 is the only free parameter to be fixed by data and here r = jri — in the potential 
picture. To convert into the momentum-space, we have 



I d Pl0x(Pl')0/(Pl)( n . _ + To— (30) 



mj — Mj^ J PI, + mv Pi+ nil V Fi'Pi 



and 



hv = / Pi </>/(?! )(-o—r 0— >\ 0—0 ' (31) 

mi - mx' J Py + my p\ + mi \j PyPi 

Ha, = ^ / d^pi(t)*x{pi')(l}i{pi) X 

mi + mx' J 



I (pI, + mv){p1 + mi) 

(p;,+miO(p? + "^i)^V ' ^ ' 

2{mi + m,x')^, iimimx* [3 Pi' 

= 3mj + mj„ - imi~mx^){3mj + mj,.) J ^'^^^P^^'^'^P'^^^^T^ + 

Pi y I iPy + mv ) {p°i + TOl) ^^^^ 



p^ + mi y PyPi 

where the 4'(x.x*) are wavefunctions of "0(1^, 2s) in momentum space, i.e. the Fourier trans- 
formed (^8|) and (^9|), in the expressions, pi and py denote the 3-momenta of the quarks which take 
part in the reaction in the initial and final mesons, while mi, my are their masses respectively, p^ 
and p° correspond to the third and the zero-th components of the concerned 4-momenta. In the 
helicity-coupling picture, all momenta of the mesons are along z, so 

PI = \Pil PxiX') = ±\PX{X')\, 

but the quark momenta can be along any directions. In the CM frame of the decaying meson 

— 

Pi = and \px{x-') \ — ' 2m^'"^*' 1^ ~ 0, thus one has 

pi+P2=pi = {M,6), and pv + P2' ^ Px{X') = {Px(X'-)^^iO^Px(X'))- 



The resultant formulae look quite different from that given in ref . Pq| , but as a matter of fact, 
as \p\ ^ M, they coincide with each other. 



Substituting all the information into eqs.(ElJ) and (22), we can have the final numerical results. 



III. The numerical results 

In the whole calculations, only /3 is a free parameter and one can fix it by the energy-minimum 
condition 

dE d<H> 
'dp " p " ■ 

Then one obtains 

where fi is the reduced mass and a is an average radius of the meson. There can be an uncertainty 
for a and fj,, it does not affect our qualitative conclusion even though indeed the numerical results 
can be declined by a few ten percents. (see below). 

Even though fo is not well measured yet, there are some reasonable estimated values, so we 
take fo = 0.15 GeV and fs = 0.125 GeV Numerically we use 

U = 0.132, fK = 0.161, fp = 0.212, fK' = 0.221 

in GeV. 

By the well-measured value as(Af|) = 0.118 ||l|, we have as{mb = 5 GeV) — 0.203, as{mc = 
1.5 GeV) = 0.265, and 

4^' = 1.26, ) ^ -0.51 

ci""^ = 1.10, cp' = -0.23. 

Our result is fully consistent with ref. |^ obtained in terms of RGE. 

It is also noted that since rric is not very large, one can expect, the real values of €{^2 may 
deviate from that predicted by the perturbative QCD calculation, for example, it is claimed that 
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a set of cf^ = 1.26 ± 0.04 and = -0.51 ± 0.05 can fit data better. However, in below we will 
rely on the perturbative QCD and use the values obtained by RGE. 
The corresponding a\^'^^ and a".^'^^ would depend on ^ of eq.(||). 

For the radially excited V'2s states, we wiU take Mo{2s) w 2.4 GeV and Mt,{2s) = 1.0 GeV, 
Mk{2s) = 1.4 GeV. By eq.(||), we fix 

13b = 0.5, (3d{Is) = 0.45, I3d{2s) = 0.39, /3^(ls) = 0.3, /3^(2s) = 0.26, %(ls) = 0.4, /3k(2s) = 0.34 



in GeV. All the parameters are obtained according to eq.(34) and (^5|). 
Numerically, we have 



D+ 




(5 = 

5 = -0.5, (36) 
(5 = -1. 



It seems that the 5 — —1 solution suits the data on D-decays better than other 5 values and this 
conclusion was also predicted by Stech ct al. a long while ago ||2^ . 

For the B-case, without considering the ■02s excited state contribution, we have 

(5 = 0, 

(5 = -0.5, (37) 

(5 = -l, 



Fro 




If one looks at (5 = — 1 which is consistent with that obtained in D-decays, the ratio is close to 
0.5 as expected (see the introduction). When we take into account the contributions from the ^2s 
excited states, the whole result is modified as 

(5 = 

(5 = -0.5 (38) 

(5= -1 



Fflo 




this result is very consistent with the data on the lifetimes of both D and B-mesons. We will 
discuss this result in next section. 



IV. Conclusion and discussion 

B and D mesons all contain a heavy quark and a light one, we have every reason to believe that 
they have similar characteristics. Indeed a symmetry between b and c quarks (B and D mesons) 
is confirmed by phenomenology. However, one obvious discrepancy that T£)± ~ 2t£)0 while 
Tg± ^ tbo implies some distinction between B and D mesons. 
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There have been ahernative ways to interpret the hfetime difference of B and D. For example, 
Bander, Silverman and Soni |^ suggested that the reaction D'^ s + d+ gluon as a source for the 
difference in the lifetimes of and and in other way, one can suppose that the factorization 
factor 6 can be different for B and D or the signs of a2 can change etc. However, if we consider 
similarities between B and D, it is natural to accept an assumption that 6 would not be too declined 
in B and D cases. In literature in D-physics, S is very close to —1 and our results confirm this 
allegation. Cheng found |^ that r2 = —0.67, —(0.9 — 1.1) for D — > Ktt, K*7t respectively, where 
our S — {Nc/2)r2, it indicates that S ^ —1. But to fit B-decay data, Cheng concluded r2 = +0.36 
which drastically deviates from the parameter for D-decays, so one would ask how it could be so? 

Instead, we accept the assumption that a symmetry between b and c holds and c[^'^\ ^2^'^^ 
can be derived with the RGE. Meanwhile we also notice that since B-mesons are much heavier than 
D-mesons, there can be radially excited states ip2s and ip2s decay products in B-decays, but not 
for D-decays. The ^2s states may cause the hadronic matrix elements to be in opposite sign to the 
ipis final states and it would result in a change to make r^i ^ Tgo . Obviously, it is determined 
by an overlapping integral between wavefunctions of the final and initial mesons. Our numerical 
results show that the integrals for 'ip2s and ipis can have opposite signs depending on the parameter 
(3. Our /3— values are reasonably determined by data, even though not very accurate. We show 
that as (5 ~ —1, as taking into account the contribution from "02^'^ as well as tp^s'^ , approximately 

Our mechanism is in parallel to the PI effects discussed by some authors Qjij. It is based 
on the common knowledge that as long as all the exclusive channels (in fact, the main ones) are 
summed up, the total width should be obtained, i.e. equivalent to the inclusive evaluation. Thus 
in our picture an interference between the decay products of the b (c) quark and the light one is 
automatically considered via the ai and 02 interference. 

Since, indeed, we only consider the most Cabibbo favorable channels to estimate the lifetimes, 
there can be contributions from the rare decays and the numerical results can deviate a bit, but 
in general the same mechanism proposed by Close and Lipkin can apply. Hence the rule is the 
same to all channels, namely always contributes as well as ipis: our results seem sufficiently 
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convincing. As a matter of fact, the 7/125 is still light enough and there is large phase space available 
for B, but on contrary, not for D- meson. 

For evaluating the hadronic matrix elements, we use the non-relativistic quark model. Even 
though the model is approximate, our qualitative conclusion docs not change. 

Surely, we can make the ratios of lifetimes for D and B mesons perfectly coincide with the data 
by carefully adjusting the [3 values in the wavcfunctions. However, since there are many uncertain 
factors such as the contributions of the rare decays, non-relativistic form of the wavcfunctions and 
the factorization factor 5 etc, which make a very accurate evaluation impossible, so only adjusting 
(3 value to fit data seems not necessary. In fact, as the most important point, one can draw 
a qualitative conclusion confidently that the the contribution of %l^2s is important to B-decays, 
namely the puzzle of the lifetimes of B and D mesons can be reasonably explained away by its 
participation. 

It is important to notice that not only the lifetimes of B- meson is in contrary to our knowledge 
based on the perturbative QCD and D-physics if the ip2s contributions is not taken into account, 
but also similar puzzles exist at many channels of B-meson decays. It is that the value of 02 is not 
universal pH and its sign is also uncertain. It is hard to understand. So we hope that by taking into 
account of the contributions, all the discrepancies may get a reasonable explanation. Because 
the relatively heavy il)2s is still light to B-meson and does not affect its phase space integration 
very much, so maybe in measurements of exclusive channels, certain ■02s with the same quantum 
numbers as gets mixed in and is not well tagged out. It causes the superficial discrepancy. 
To carefully and thoroughly investigate the influence and effects of possible ■023 decay products in 
B-decays is the goal of our next works. 
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Figure Captions 

Fig.l(a)-(d). The quark diagrams for the non-leptonic decays of B and D mesons (here we take 
D — > Kit as an example). 
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